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SUMMARY 


The elastic constants of a fiberglass epoxy unidirectional 
composite are determined by measuring the phase velocities of 
longitudinal and shear stress waves via the through transmission 
ultrasonic technique . The waves introduced into the composite 
specimens were generated by piezoceramic transducers. Geometric 
lengths and the times required to travel those lengths were used 
to calculate the phase velocities. The model of the transversely 
isotropic medium was adopted to relate the velocities and 
elastic constants. 
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INTRODUCTION 


Filamentary composites have proved their usefulness in 
modern aircraft, spacecraft, land and water vehicles, and many 
other structural applications. 

One of the special characteristics of these materials is 
their anisotropic elastic behavior, that is, their elastic 
properties which change with orientation. Knowledge of such 
properties is essential both for optimal design and for 
evaluation of performance. Elastic properties are usually 
represented by a set of constants which can be determined by 
conventional stress-strain measurements in different loading 
configurations (as, for instance, in standard tensile tests) or 
by ultrasonic methods. 

Ultrasonic testing is a nondestructive evaluation (NDE) 
technique in which the propagation characteristics of stress 
waves in a material are determined . The wave transmission 
characteristics may then be related to such factors as the 
elastic properties of the material or the presence of flaws. 
Ultrasonic methods have been extensively used for the 
determination of velocities of stress waves, elastic constants of 
the stiffness matrix, fiber fraction, degree of fiber 
orientation, presence of voids, and attenuation and damping 
parameters [1—6] . One of the main advantages of ultrasonic 
methods, as they are essentially nondestructive, is their 
suitability for the evaluation of !, in service" components . 

The purpose of this work is to establish the experimental 
values of the velocities of ultrasonic waves of a specific 
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composite laminate of fiberglass epoxy in which the fibers are 


continuous 
the values 
determined . 
purposes 
performance 


and unidirectional. From the measured velocities, 
of elastic constants of the stiffness matrix are 
Such parameters can then be used both for design 
or as standards for the NDE evaluation of the 
of components. 
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TRANSVERSELY ISOTROPIC COMPOSITE 

In filamentary composite materials, the presence of a 
distinctive fiber orientation gives rise to elastic anisotropy. 
Unidirectional fiber composites consist of a parallel set of 
cylindrical fibers embedded in an otherwise homogeneous matrix 
material. If the constituent materials are isotropic, the 
resulting composite will be classified as transversely isotropic, 
or hexagonally symmetric. 

A transversely isotropic material can be formed from 
prepreg sheets of continuous parallel fibers that are layered to 
form a unidirectional composite [7]* A schematic of this 
configuration is shown in Fig.1. When the transversely isotropic 
model is adopted, the fibers are assumed to be randomly arranged 
in the xy plane. 

To characterize the mechanical behavior of transversely 
isotropic materials, five independent elastic constants are 
required [8] . 
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BASIC PRINCIPLES 




The elastic properties of a filamentary transversely 
isotropic composite are evaluated by measuring the phase 
velocities of longitudinal and shear plane waves propagating 
through the material at various angles with respect to the 
fiber direction. 

Depending on the direction cosines chosen for the normal of 
the propagating wave, the equations yielding the five elastic 
constants in terms of velocities can be obtained. These equations 
have been derived for transversely isotropic materials [9] . For 
a cartesian coordinate system x, y and z in which z coincides 
with the fiber direction in the transversely isotropic material, 
the relations between the elastic constants of the stiffness 
matrix and phase velocities are summarized in Table 1 . The 
propagation of waves in any direction or polarization direction 
of a wave is represented by its direction cosines with respect to 
the x, y, and z directions (for example (1,0,0) defines the x 
direction). P denotes an imposed longitudinal wave, for which 
propagation and polarization directions are coincident; SH 
denotes an imposed shear wave with polarization in the y axis ; 
and SV denotes imposed shear wave with polarizaton in the x-z 
plane . 


The constants C „, C „ . 

11 33 44 


and 0^^ can be determined by 


direct measurement of the velocities of appropriate waves 
propagating in the principal directions x,y and z. The 


remaining constant, C 
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can be determined 


for 


wave 


propagating along an intermediate direction between the axes x 
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and z . (Table 1 shows the equations for some intermediate 

directions, defined by their angles with respect to the z axis.) 

The constant C.„ can be calculated from C„. and as 

12 1 1 66 

as C 12 = C 11~ 2C 66 * 

The model adopted describes the composite as homogeneous, so 
the density is assumed to be constant. The model is valid for 
stress waves having wavelengths that are large compared with the 
thickness of the prepreg sheets. Previous studies (10,11 ] have 
shown that waves can be scattered by reflection if the 

wavelengths have dimensions comparable to the thickness of the 

lamina or the fiber diameter. Effects of dispersion may be 
present in such cases, and a better model ought be a layered 
medium model instead of the homogeneous model. 

For the homogeneous model, the correspondence between the 
elastic constants of the stiffness matrix and the magnitudes of 
extensional moduli, shear moduli, and Poisson’s ratios are shown 
in Table 2. 

For composites with symmetries other than hexagonal, the 
procedure of measuring velocities can also be used, but in these 
cases a different set and number of independent elastic constants 
must be determined . It would be possible to determine the 
symmetry axes (if any ) of any specimen by making sufficient 
measurements. This would involve a large amount of experimental 
work and analysis which could be greatly reduced if some prior 
information is known about the symmetries of the material. 


5 



EXPERIMENTAL PROCEDURES 


Procedures and Equipment 

The measurements were made by the through transmission 
technique using specimens of different lengths for each of the 
measurement directions. A schematic of the measuring system is 
shown in Fig. 2 . The system consisted of a pulsed oscillator 
(Arenberg model PG-652) for generating sinusoidal waves , a low 
frequency inductor (Aremberg model LFT-500 ),two broadband (0.1 
to 3.0 MHz) transducers (AET model FC-500) for longitudinal waves 
with approximately flat sensitivity of -85 dB (relative to IV/ 
Bar), two transducers (Panametrics model V154) for shear waves 
with nominal resonance of 2.25 MHz , an oscilloscope (Tektronix 
model 455) > auxiliary attenuators , and a pneumatic fixture 
that provides constant pressure between transducers and specimen 
interfaces. Couplants AET SC-6 and Panametrics SWC were used at 
interface of transducers and specimens for longitudinal and 
shear tests respectively. 

The transmitting transducer is excited with a tone burst 
voltage signal of amplitude 100 Volts (peak to peak); this 
signal imposes a stress wave (assumed plane) onto the specimen. 
After travelling through the specimen, the signal is captured by 
the receiving transducer and displayed on the oscilloscope. 

The velocity measurements of ultrasonic waves in the 
through transmission configuration may be affected by many 
parameters, including the transducer-specimen interface pressure, 
the mechanical properties and thickness of the couplant, the 
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the 


alignment between the specimen and the transducer, 
nonparallelism of the faces of the specimen and the inherent time 
delay of the electrical system. To avoid these effects, the 
following steps were taken: 

a) The specimen was very carefully centered on the transducer 
face . 

b) The interface pressure was maintained at 4.2 Bar during 
all the tests. This value is slightly above the saturation 
pressure [5] and guarantees a stable output. 

c) The specimens were machined within established tolerance 
levels for parallelism between faces . 

d) To evaluate the velocity , the time shift between the same 
individual cycle of the input and ouput signals was measured. A 
typical example of the input and ouput signals is shown in Fig.3« 
A few leading cycles were not used because of their transient 
nature. In order to cancel the inherent time delay due to the 
electrical system, otherwise identical specimens of different 
lengths were used. If the time shifts for specimens of lengths 
L. are t^, respectively, the phase velocity v can be expressed as 

v = (L 2 - L 1 )/(t 2 - t.,) (1) 

The tests were conducted at different input frequencies: 
0.28, 0.50 ,1.0, 1.5 and 2.0 MHz . The selected frequency 

range depended on the frequency range of operation of the 
transducers and on the purpose of application of the values 
measured . 
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Specimens 


The material used was 3M Scotchply 1002, continuous glass 
fiber (type E ) in a 165°C (329°F) curing epoxy matrix. The 
composite was cured in a heated press at 0.69 MN/m (100 psi) 
according 3M specifications [12] . The resulting resin content 
was 36 percent by weight. Initially a 25.4 mm (1 in) thick plate 
with dimensions 355.6mm x 279*4 mm (I4in x12in) was fabricated. 
After removing the edges (approximately 10 mm ), the specimens 
were cut from the plate with a diamond saw. The positions of the 
specimens with respect to the longitudinal axis of the plate (z) 
are shown in Fig .4 . 

The cross-sectional dimensions were, based on previous 
experience [5,13] » smaller than the diameter of the transducers 
faces. All cross-sections were square with 12.7 mm (0.5 in) 
sides . 

The length to be traversed by the waves, ranged from 
5mm (0.19in) to 29.6 mm (1.06in). These dimensions were limited 
by the geometry of the original plate and the positions of the 
specimens with respect to the longitudinal axis of the plate . 
Attention was given to the degree of parallelism of the final 
surfaces rather than to the dimensions , since the last could be 
measured a posteriori for the velocity calculations. The degree 
of parallelism (along the length) was maintained better than 
+.013 mm (+.0005in) on all specimens. 

Since the velocities are to be measured in 5 different 
directions with respect to the longitudinal axis of the plate 
(0°, 30°, 45°, 60°, 90°), at least 10 specimens (2 specimens 
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per direction) would be required to make all the measurements . 
However , to allow for possible errors and variations in the 
material , three specimens of different lengths were used for each 
of the directions. 


9 


RESULTS 


The times measured for the waves to travel through a certain 
thickness of material were found to be independent of the 
frequency for each of the investigated directions. This finding 
shows that the material is nondispersive in the frequency range 
used. The wave speeds were calculated from the time measurements 
in each direction, and the results were averaged. The confidence 
level of time measurements is 95 percent. The results are shown 
in Table 3« 

The largest values for the velocities occurred for the 
imposed longitudinal waves, P , followed by the shear waves with 
polarization in the xz plane, SV, and finally by the shear waves 
with polarization in the y direction, SH. 

For the purpose of comparison, the values of velocities 
along the principal directions were also calculated from layered 
media theory [14] and are also shown in Table 3« Table 3 shows 
also the percentage differences between the experimental and 
layered media theory values. 

From the experimental values of the velocities , the elastic 
constants of the stiffness matrix were calculated using the 
expressions in Table 1 . The value used for the density was 
1850 kg/m^ according to manufacturer's data [12]. The final 
results for the elastic constants of the stiffness matrix are 
shown in Table 4 with the corresponding elastic moduli and 
Poisson's ratios shown in Table 5. Observe that the 
measurements along nonprincipal directions of propagation ( such 
as the directions 30°, 45° and 60° with respect to z ) are 
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supposed to give values of velocities that can be used for the 

calculation of the same constant (see Table 1). The 

apparently excessive number of measurements was done because for 

directions other than the principal ones (x, y and z) , the error 

in the time measurements can be large. This is due to the fact 

that the directions of plane waves originally imposed onto the 

material are modified by a deviation effect [15] , possibly 

introducing errors in measured times. The final value of the 

constant was obtained by averaging the individual values 

of C , calculated for each of the nonprincipal directions 
13 

(30°, 45°, 60°) measured. 

In general, whenever the same elastic constant of the 
stiffness matrix could be calculated from measurements along 
different directions, the final value was obtained by averaging 
the values obtained for that constant in each individual 
direction . 
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CONCLUSIONS 


The measured phase velocities indicate that the fiberglass 
epoxy material tested is nondispersive for frequencies between 
0.28 and 2.0 MHz. 

The comparison of the experimental results with the results 
of layered media theory show reasonable agreement for the values 
of velocities of longitudinal waves. For shear waves, the 
percentage differences between some of the experimental and 
theoretical values is large. 

The values of the velocities were used to calculate the 

elastic constants of the stiffness matrix, elastic moduli and the 

Poisson's ratios for the fiberglass epoxy. The values found are 

elastically compatible. For the compatibility of elastic 

properties the Poisson's ratios v and v must obey the 

xy zx 

following condition [16] 

-1 * v £ 1 

( 2 ) 

V *x * <W 1/2 

where E and E are the extensional moduli of the composite in 
x z 

the x and z directions, respectively. 
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TABLE 1 Relations Between Elastic Constants of Stiffness Matrix and Phase Velocities. 
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represents the density and v the phase * veloci ty in the propagation direction 


TABLE 2 Relations between Elastic Constants, Elastic 
Moduli and Poisson's Ratios for Transversely 
Isotropic Medium. 


E = 

X 

< C 11- C 12 

^ C 11 C 33 + C 12 C 33 

-2C 13 2 ) / (C 1lC „ - C 13 2 ) 

E = 

y 

E x 



E = 
z 

(C 11 C 33 + 

C 12 C 33 ” 2C 13 ^ ^ 

^ C 11 + °1 2 ^ 

II 

M 

(C 12 C 33- 

C 13 ^ ^ ^ C 11 C 33~ C 

1 3 2 ) (*) 

v yx“ 

v xy 



v xz* 

C 13 1 < 

C 11 * C 12 > 


v zy 

v zx 



v zx 

* - C (C 

13 v 

12 * °11^ ^ ^ C 11 

C 33 “ °13 ) 

v zy 

s v 2X 



G 

xz 

* C 44 



G 

yz 

= G xz 



G 

xy 

* (c 11 - c 

12 >/ 2 



(*) Characterizes the contraction in the x direction when tension 
is applied in the y direction. 
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TABLE 3 Velocities for Unidirectional Fiberglass Epoxy 


1 

Direction | Imposed 
of | Wave 

Propagation, 1 Type 

i 

i 

i i 

Velocity (m/s) 

Experimental j Layered 

| Media Theory 

i 

j Percentage j 
j Difference | 
i Experimental | 

i 

z 

i P 

i 

4692.8 

1 

1 

4344.7 

j 7.4 j 

30° (*) 

i p 

4181 .0 

1 

1 

1 

— 

! — ! 

45° 

i p 
1 

3558.7 

1 

1 

— 

i — | 

60° 

! p 

2940.9 

1 

1 

— 

i — i 

X 

i P 

2391 .9 

1 

1 

2300.8 

i 3.8 j 

z 

j SH 

i 

1550.4 

1 

1 

1380.6 

i 10.96 i 

30° 

| SH 

1521 .1 

1 

1 

— 

! — ! 

45° 

j SH 

1480.6 

1 

1 

1 

— 

! — ! 

60° 

j SH 

1399.5 

1 

.1 

— 

! — j 

X 

SH 

1343.0 

i 

1 

1380.6 

i -2.7 ! 

z 

i SV 

1586.3 

1 

1 

2835.0 

! -78.72 | 

30° 

| SV 

1704-7 

1 

1 

— 

1 — ; 

45° 

! SV 

1836.1 

1 

— 

1 — ! 

60° 

| SV 

1904.2 

1 

1 

— 

j — | 

i! x 

| SV 

1506. 3 

1 

1 

1380.6 

1 

1 

1 

» 

1 00 

1 • 

1 01 

1 ^ 
1 
1 
1 
1 


(*) All angles are with respect to z axis. 
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TABLE 4 Values of Elastic Constants of the Stiffness Matrix 


I Elastic Constant 

! Experimental Value (GN/m ) 

i 


c n 

-X 

O 

• 

Ul 

03 

i i 


c i2 

! 4-098 

t 

i 


j C 13 

1 4-679 

i 

i 


1 c ” 

j 40.741 

i 

i 

i 


| C 44 

i 4-422 

i 

i 


| C 66 
1 

! 3-243 

i 

i 
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TABLE 5 Values of Elastic Moduli and Poisson's Ratios 


i ! __ ! 

| Elastic Moduli \ 

land Poisson ratios! 

■ i 

Ezperimental Values 


j Eztensional 

E i 

z 1 

i 

i 

37.759 

GN/m 2 


j Eztensional 

E z | 
1 
1 

8.785 

GN/m 2 


| Poisson 

V xy | 
1 

0.354 



j Poisson 

V i 

XZ ! 
1 
1 

0.319 



| Poisson 

V 

Z X ! 
1 
1 

0.074 


} 

| Shear 

G ! 
xy i 

i 

3.243 

GN/m 2 


j Shear 

g i 

XZ J 
I 
1 

4-422 

GN/m 2 
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Schematic of experimental system for through transmission 
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